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Abstract. In this paper we describe some analytical properties of the R}/ m law proposed 
by Sersic (1968) to categorize the photometric profiles of elliptical galaxies. In particular, 
we present the full asymptotic expansion for the dimensionless scale factor b(m) that is 
introduced when referring the profile to the standard effective radius. Surprisingly, our 
asymptotic analysis turns out to be useful even for values of m as low as unity, thus 
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' providing a unified analytical tool for observational and theoretical investigations based 
on the R}l m law for the entire range of interesting photometric profiles, from spiral to 
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elliptical galaxies. 
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1. Introduction 

After its introduction as a generalization of the i? 1 / 4 law (de Vaucouleurs 1948), the so- 
called Sersic law (Sersic 1968) has found a variety of applications. On the observational 
side, it has been used as a tool to quantify the non-homology of elliptical galaxies (see, 
e.g., Davies et al. 1988; Capaccioli 1989, hereafter C89; Caon, Capaccioli & D'Onofrio 
1993; Young & Currie 1994; D'Onofrio, Capaccioli & Caon 1994; Prugniel & Simien 
1997, hereafter PS97; Wadadekar, Robbason & Kembhavi 1999). In addition, it has been 
applied to the description of the the surface brightness profiles of galaxy bulges (see 
Andredakis, Peletier & Balcells 1995; Courteau, De Jong & Broeils 1996). One research 
area where the usefulness of the Sersic law as a statistically convenient description has 
been exploited is that of the Fundamental Plane of elliptical galaxies (Graham et al. 1996; 
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Ciotti, Lanzoni & Renzini 1996; Graham & Colless 1997; Ciotti & Lanzoni 1997; Graham 
1998). On the theoretical side, it has been the focus of several general investigations (see, 
e.g., Makino, Akiyama & Sugimoto 1990; Ciotti 1991, hereafter C91; Gerbal et al. 1997; 
Andredakis 1998). 

According to such law, the surface brightness profile is given by 

I(R) = I Q e- b n 1/m , (1) 

where n = R/ R c , m is a positive real number, and 6 a dimensionless constant such that 
R c is the effective radius, i.e., the projected radius encircling half of the total luminosity 
associated with I(R). In a broad statistical sense, it is found that bright ellipticals are 
well fitted by the Sersic law with m around 4, dwarf ellipticals and galaxy disks with m 
around 1, and finally bulges and intermediate luminosity ellipticals with 1 < m < 4. For 
some galaxies, a value of m even higher than 10 has been found (e.g., see NGC 4552, 
Caon et al. 1993). 

The projected luminosity inside the projected radius R is given by 

L(R) = 2tt j* I{R')R'dR' = IvRl^^mM 11 ™), (2) 
where (for a > 0) 

l{a,x) = f e-H a - x dt (3) 



Jo 

is the (left) incomplete gamma function. The total luminosity is then given by 

L = I ^^T(2m), (4) 

where T(a) = 7(0, 00) is the complete gamma function. From the definition of R e it 
follows that 6(m) is the solution of the following equation: 

7 (2m, & ) = M. (5) 



2. Asymptotic expansion 

Unfortunately, Eq. (5) cannot be solved in explicit, closed form, and so it is usually 
solved numerically^. This is inconvenient for a number of observational and theoretical 
applications. The exact values of b(m) are recorded in Table 1 for 1 < m < 10. For the 
de Vaucouleurs law, m = 4 and 6(4) « 7.66924944. Interpolation formulae for b(m) have 
been given in the literature, namely b ~ 1.9992m— 0.3271 by C89 (as reported by Graham 
& Colless 1997), b ~ 2m - 0.324 by C91, b ~ 2m - 1/3 (for m integer) by Moriondo, 
Giovanardi & Hunt (1998), and the "numerical solution" b(m) ~ 2m— 1/3 + 0.009876/m 
For m = 1, i.e., the exponential profile, the solution can be formally expressed using the 
Lambert W function, as 6(1) = -1 - W(-l, -l/2e) = 1.678346990... 
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by PS97. These expressions provide an accurate fit in the range 0.5 < m < 10; curiously, 
their leading term is linear in m, with a slope very close to 2. In the following we show 
that this behavior results from a general property of the gamma function. 

Prompted by Eq. (5) we now address the following: 
Problem Solve for x 

7 ( tt>a r) = -Li, (6) 
for given a > 0. 

Because no explicit solution in closed form is available, we will focus on the asymptotic 
expansion of x{a) for a >> 1. In fact, it is well known that in many cases asymptotic 
expansions turn out to give excellent approximations of the true function even for rel- 
atively small values of the expansion parameter. The starting point of our study is the 
asymptotic relation (see Abramowitz & Stegun 1965) 



V a 



1 1 139 571 163879 

Yhx + 288a 2 ~ 51840a 3 ~ 2488320a 4 + 209018880a 5 



This is the Stirling formula^, which is known to be associated with a relative error smaller 
than 3 x 10 -6 already for a = 2. 

Let us now introduce the sequence 

n-l 

X n = CI + V — r (8) 

i — < a 

k=0 

with so = a, so that x n +\ — x n + c n /a n . Here Ck are coefficients (to be determined 

at a later stage), independent of a. Then we start by proving the following asymptotic 

results, applicable for a >> 1. 

Lemma 1 The following asymptotic relation 

w \ 00 p(0) 

iM-^ + e-a"^^, (9) 

fc=0 

holds, where Pjf are rational numbers. 

The validity of Eq. (9) can be established by means of a standard asymptotic expansion 
(e.g., see Bender & Orszag 1978, Bleinstein & Handelsman 1986) of the integral 

*»°>-«-«- °/; ^- + + r (l+ ^ 

In fact, the argument of the integral is the same as that of the integral representation of 
r(a). In both cases the stationary point for the exponent occurs at s = 0, but for T(a) 
the stationary point is in the middle of the domain of integration, because the integral 



2 The derivation of this formula can be found in standard textbooks. The coefficients appearing 
in the asymptotic expansion of lnF(a) for a — * oo can be expressed in terms of the so-called 
Bernoulli numbers; see, e.g., Arfken & Weber 1995, Chapts. 5 and 10. 
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extends to oo (instead, for the integral in Eq. [10] the upper limit is precisely s = 0). 
Thus, when we consider the power series expansion (in s) of the argument of the integral 
around the stationary point, for 7(0:, a) the even powers of s contribute exactly one 
half of their contribution to T(a), while the odd powers determine the terms in Eq. (9) 
associated with the coefficients P^ (in contrast, the odd powers do not contribute to 
r(a), by symmetry). The calculation of Pjf^ is tedious, but straightforward. 

Note that there is a "shift" of powers, by a 1 / 2 , between the two terms on the right 
hand side of Eq. (9). In particular, the second term is smaller by a factor 0(a -1 / 2 ). This 
already shows that x = a is a first approximate solution to the problem set by Eq. (6). 
Lemma 2 The following asymptotic relation 

7(a, x n+1 ) ~ -f(a, x n ) + e~ a a a f(a) (11) 
holds, with f(a) = 0(a~™~ 1 ). To leading order, f(a) <~ c n /a n+1 . 

This result easily follows from the definitions of the quantities involved (Eqs. [3] and [8]), 
which give 

rc„/a n 

-y(a,x n+1 ) =-y(a,x n ) + e- Xn / e~\t + x^^dt. (12) 

Jo 

At this point we can proceed to prove the following theorem: 
Theorem For large (real) values of a, the full asymptotic expansion of the solution to 
the problem posed by Eq. (6) can be expressed as 



x(a)=a + J2^, ( 13 ) 

71=0 

where 

Cn = -P^\ (14) 

and the coefficients can be calculated by iteration on the relation 

00 (n) 

T (a,x„)~rH/2| e -V^-^. (15) 

k—n 

The proof is obtained by induction. In fact, Eq. (9) shows that the statement is true 
for n = 0, with the coefficients P^ available from the asymptotic analysis outlined in 
the proof of Lemma 1. If we now refer to the result of Lemma 2, with the leading order 
expression for f(a), and assume the statement (related to Eq.[15]) to hold true for x n , 
we find 

_k—r, 



j{a,x n+1 ) 



p(™) 
1 1 



(16) 



In other words, the statement is found to hold true also for x„ +1 , provided c„ = — Pn n \ 
as required by Eq. (14). The method thus provides a way to systematically improve the 
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approximation to x(a) by means of n steps, leading to an estimate x n ; step by step the 
possible presence of undesidered "shifted" (odd) terms is eliminated and the process leads 
to the complete determination of the coefficients defining the asymptotic series (13). At 
a given level n of desidered accuracy the coefficients depend on the values P^ for 
i = 0, ...,n — 1. 

The explicit computation yields: 

_1 _8_ 184 1048 17557576 _ 5 

X[a) ~ a 3 + 405a + 25515a 2 + 1148175a 3 15345358875a 4 + [a >' [ ' 

The first two terms can be easily checked using standard general formulae for the leading 
terms of the relevant steepest descent asymptotic expansion. 



3. Analytical properties of the Sersic law 

Therefore, the first terms of the asymptotic expansion of b(m) (for real m) are 

j . . 1 4 46 131 2194697 n , _ 5 . „„. 

&(m) ~ 2m ~ 3 + i05^ + 25515^ + 1148175^ ~ 30690717750m 4 + ° (m " } ' (18) 

Equation (18) now clearly explains the value of the interpolation formulae found earlier 
(C89, C91, PS97); note that 4/405 = 0.0098765.... How many terms in the asymptotic 
expansion are required to obtain a better representation of b(m) when compared to the 
previously introduced interpolations? 

We have computed the relative errors of the various expressions with respect to the 
true value of b(m) (obtained by solving numerically Eq. [5] with a precision of 20 sig- 
nificant digits) for integer values of m in the range 1 < m < 10, and the results are 
reported in Table 1. The first result is that using the first four terms of the expansion 
the true value of b(m) is recovered with a relative error of 6 x 10~ 7 for m = 1, and 
4 x 10~ 9 for m = 10, i.e., the asymptotic expansion so truncated performs much better 
than the formulae cited previously. Obviously for larger values of m the error becomes 
correspondingly smaller. The second somewhat surprising result is the fact that Eq. (18) 
is already very accurate for m as small as 1. This allows us to include, within the reach 
of the present analysis, the case of exponential profiles. A third point that we have noted 
is that, for fixed m, there is an optimal truncation of the asymptotic expansion, beyond 
which, as is well known in the general context of asymptotic analysis, increasing the 
number of terms in the expansion does not improve the accuracy of the estimate. For 
example, for m = 1, the optimal truncation occurs at the fourth term, for which the 
attained relative error is 6 x 10~ 7 . For simplicity in the following part of this Section 
we will record a number of interesting analytical expressions restricted to their leading 
order. Of course, the asymptotic analysis provided here would allow us to give explicitly 
any higher order term, not shown below, if so desired. 
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Table 1. True values and relative errors on b(m) for integer values of to, using the C89, 
C91, and PS97 formulae. As(4) is the relative error for the asymptotic expansion given 
in Eq. (18), truncated to the first four terms. 



m 


b(m) 


C89 




C91 




PS97 




As(4) 


1 


1.67834699 


4 x 10" 


■3 


10~ 3 




10~ 3 




6 x 10" 7 


2 


3.67206075 


2 x 10" 


-4 


10 -3 




1 1 1 — 4 

10 




10 -6 


3 


5.67016119 


6 x 10" 


-5 


10~ 3 




4 x 10" 


-5 


4 x 10 -7 


4 


7.66924944 


6 x 10" 


-5 


9 x 10" 


-4 


10 -5 




io- 7 


5 


9.66871461 


2 x 10" 


-5 


8 x 10" 


-4 


8 x 10" 


-6 


5 x 10~ 8 


6 


11.6683632 


2 x 10" 


-5 


7 x 10" 


-4 


4 x 10" 


-6 


3 x 10" 8 


7 


13.6681146 


6 x 10" 


-5 


6 x 10" 


-4 


3 x 10" 


-6 


2 x 10~ 8 


8 


15.6679295 


9 x 10" 


-5 


5 x 10" 


-4 


2 x 10" 


-6 


9 x 10" 9 


9 


17.6677864 


io~ 4 




5 x 10" 


-4 


10 -6 




6 x 10~ 9 


10 


19.6676724 


io~ 4 




4 x 10" 


-4 


9 x 10" 


-7 


4 x 10~ 9 



3.1. Total luminosity and central potential for a spherical system with an R 1 /™ 1 
projected luminosity profile 

From Eqs. (4) and (18) the total luminosity is found to be 

L = I Rl^-T{2ra) ~ hR 2 ^' 2 e x ' z e~ 2m ^[1 + O(to^)], (19) 

where we have used the fact that b 2m - e- 1 / 3 (2m) 2m [l + 0(m- 1 )]. Following C91, the 
central potential of the spherically symmetric density distribution associated with the 
Sersic law is given by: 

*o = -4^ffl^l ,20, 

Here M is the total mass of the system, and the mass-to-light ratio is taken to be 
constant. Thus, an asymptotic estimate at given I , R e is 

$ ~ j oHb 2»/V/V/ 8 ^[1 + 0(m-% (21) 

We note that, for to = 1 and m = 4, a truncation of Eq. (19) to the leading term is 
characterized by a relative error of 5.7 per cent and of 1.5 per cent, respectively. The 
corresponding truncation on Eq. (21) is associated with a relative error of 8.6 per cent 
and 2.3 per cent. 

The proper normalization of the R 1 /™ profile, to be considered for a case with given 
scales L and R c , is 

L h 2m e - b v 1/m 
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which thus provides the useful quantity 

I. = I{v = D-^^=[l + 0(m-% (23) 
so that 

HM 9 3 / 2 /p\ m 

For to — 1 and to = 4, a truncation of Eq. (23) to the leading term is characterized 
by a relative error of 7.3 per cent and of 1.8 per cent, respectively. The corresponding 
truncation on Eq. (24) is associated with a relative error of 3.1 per cent and 0.8 per cent. 

4. Conclusions 

In this paper, the full asymptotic expansion for the dimensionless scale factor b(m) ap- 
pearing in the Sersic profile has been constructed. It is shown that this expansion, even 
when truncated to the first four terms as 
14 46 

Urn) = 2to 1 1 - (25) 

V ; 3 405m 25515to 2 1 ; 

performs much better than the formulae given by C89, C91 and PS97, even for to values 
as low as unity, with relative errors smaller than ~ 10~ 6 . The use of this simple formula 
is thus recommended both in theoretical and observational investigations based on the 
Sersic law. With the aid of this formula, we have been able to clarify a number of 
interesting properties associated with the Sersic profile. The additional material presented 
in Appendix A can be compared to the simple power law i?~ 2 , often used in the past to 
fit the photometric profiles of elliptical galaxies. 

Acknowledgements. We would like to thank F. Simien for several useful suggestions. This work 
was supported by MURST, contract CoFin98. L.C. was also partially supported by ASI, contract 
ASI-ARS-96-70. 

Appendix A: Remarks on a power law approximation 

The surface brightness profile given in Eq. (1) is sometimes expressed as (see Eqs. [22]- 
[23]) 

I(R) ^ I c e b(1 -^ 1/m K (A.l) 

In this case, for a fixed location r\ > 0, we may consider an asymptotic expansion of the 
surface brightness profile at constant I e for to >> 1. This is obtained by introducing a 
stretched radial coordinate £ = In r\ and by noting that 

6(1 - V 1/m ) = 6(1 - e« /m ) = -2t(l--L + —|— .\(i + -L + J—..} (A.2) 
v ' ; v ' s V 6to 405m 2 \ 2m 6m 2 J v ' 
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Fig. A.l. Ratio between the R}l m profile (normalized at 7 e ) and a properly scaled R 2 
profile, for m = 1, 4, 10, 20. 

Thus, we find 

J(i2)~ krll + Oim- 1 )}. (A.3) 

We may recall here that the photometric profiles of elliptical galaxies have often been 
described in the past in terms of power laws (see, e.g., Hubble 1930). A naive inspection of 
the first term omitted in the expansion (A2) suggests that Eq. (A3) is adequate provided 
1/3) /m | << 1. Note that the term involved vanishes at £ = and at £ = 1/3. This 
is an indication that the quality of the power law approximation is asymmetric, with a 
modest bias to the outer region. 

Consider the function /(£) = 6(1 — e^/" 1 ) + 2£. The quantity exp(/) gives the ratio 
between the i? 1 /™ profile and its R~ 2 approximation. The function / diverges to — oo 
both for £ — > — oo (i.e., r\ — > 0) and for £ — > +oo (i.e., 77 — > +00), and has a single 
maximum /m at £m, denned by the relation e^ M l m = 2m/ b, where 

/m = 6 - 2m + 2m In ~ ^- + 0(m- 2 ). (A.4) 
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Note that 2m/b is close to unity (in fact, £m ~ 1/6), i- c -, that ^u/m is small. There- 
fore, the power law profile is slightly underluminous with respect to the R 1 /™ profile 
in the radial range between the effective radius and an outer radius £ r ~ 2£m <~ 1/3, 
where /(£ r ) = 0, which coincides with the outer location identified by a previous naive 
inspection (see comment after Eq. [A3]). Outside such radial range the power law profile 
is brighter than the i? 1 /" 1 profile. In such region, where / < 0, we may ask how far (in 
radial range) Eq. (A3) applies, by studying the condition 1 — < e, i.e., |/| < e. We 
expand / for negative values of £ around £ = 0, and around £ r for £ > £ r . The range of 
applicability of Eq. (A3) is thus constrained by the condition (1 — 3e) m <?7<(l + 3e) m e 1 / 3 . 
This situation is illustrated in Fig. 1. 
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